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We extend our previous investigations of fermionic condensation in broad Feshbach resonances 
by using the two-channel model developed for narrow Feshbach resonances. We investigate two 
crossovers: the BCS-BEC crossover by changing the s-wave scattering length and the crossover from 
a narrow to a broad resonance by changing the atom-molecule coupling. At zero temperature we 
analyze, as a function of both atom-molecule coupling and s-wave scattering length, the chemical 
potential, the energy gap, and the condensate fraction of atoms. In particular, we predict the 
contribution of Cooper pairs and preformed molecules to the total condensate density along the two 
crossovers. 
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I. INTRODUCTION 

Nowadays there is a remarkable experimental control 
over ultracold alkali atoms interacting through a Fes- 
hbach resonance. Manipulation of the binding energy 
through external magnetic fields enables experimental- 
ists to evolve clouds of two-component fermionic atoms 
from the weakly coupled BCS-like behavior of Cooper 
pairs to the strongly coupled Bose-Einstein condensation 
(BEC) of molecules Q. The transition is characterized 
by a crossover in which, most simply, the s-wave scatter- 
ing length a s diverges as it changes sign ■ Recently, a 
considerable theoretical effort |4l4l2l| has been expended 
on studying the condensate fraction of Cooper pairs of 
such a tunable supcrfiuid. The behavior of tunable gases 
can differ according as the Feshbach resonance is broad 
or narrow, the former essentially describinga one-channel 
system, the latter a two-channel system [18l - |2ll ]. In two 
experiments [H|,[ll| with a broad Feshbach resonance the 
condensate fraction of Cooper pairs [l5| has been studied 
in two hyperfine component Fermi vapours of 6 Li atoms. 
The experimental data of the condensate fraction, which 
is directly related to the off-diagonal- long-range order of 
the two-body density matrix of fermions [l6l \vl\ , are in 
quite good agreement with broad-resonance mean-field 
theoretical predictions (J, [B[ at zero temperature, while 
at finite temperature beyond-mean-field corrections are 
needed @. 

In this paper we extend our previous investigations of 
fermionic condensation in broad Feshbach resonances d, 
0-H, [H| by analyzing the condensate fraction of Cooper 
pairs and the condensate fraction of preformed molecules 
in the case of a narrow Feshbach resonance. We find that 
by increasing the resonant atom-molecule coupling the 
condensate fraction of preformed molecules is strongly 
reduced in the BEC regime while the condensate fraction 
of Cooper-paired atoms grows reaching a maximum value 
close to the unitarity limit (infinite scattering length). 
Finally, for very large values of the resonant coupling we 
recover the broad-resonance regime where there are no 



more preformed molecules and the BCS-BEC crossover 
is enterely due to Cooper-paired atoms. 



II. EFFECTIVE ACTION AND SADDLE-POINT 
APPROXIMATION 

We adopt the path integral formalism [23| and consider 
a mixture of fermionic atoms and molecular bosons, in 
which the fermions, described by the complex Grassmann 
fields - 0(r(r, t), with spin a = (t,i), can be bound into 
a molecular boson described by the complex scalar field 
4>(r, t) through a Feshbach resonance (Tg. |20|. The par- 
tition function Z of the system at temperature T can be 
written as 



Z 



Vip a T>ip*-D(l>D<f>* exp (-S/K) , (1) 



where the Euclidean action functional S is given by 



with 
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dr / cFr (j*f F + Se B + &fb) 
Jv 



dr 2m 
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the Lagrangian density of free fermionic atoms, 
Jz?b = 



h— - — V 2 -2 M + e 
or 4m 



(2) 



(3) 



(4) 



the Lagrangian density of free bosonic molecules, and 

&fb = 9 (0* ^ Vt + <t> i>t) ■ (5) 

the Lagrangian density of fermion-molecule coupling. 
Notice that (5 = l/(fcsT) with ks the Boltzmann con- 
stant and V is the volume of the system. Moreover, 
the bound molecular bosons of Feshbach resonance have 
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twice the mass of the fermions and a tunable binding 
energy e . 

The action functional S is quadratic in the fermionic 
fields ip a ( r , T )i which can be then integrated out exactly 
obtaining 

Z = j V<f)Vcj)* exp (Se/h) , (6) 
I 



where S e is the effective action, given by 

rW r 

S e = -HTr[ki(&- 1 ))+ / dr / d 3 r£? B (7) 
Jo Jv 

with Sf _1 (r, r',r, r') the inverse Grccn-Nambu function, 
defined as 



"St ^ 2m v ^ / 



Within the saddle-point approximation [23j we con- 
sider a space and time-independent molecular field, i.e. 
<j>(r, t) = 0o- Without loss of generality we suppose that 
<j>o is a real number. In this way the corresponding saddle- 
point effective action function S sp obtained from Eq. ([7]) 
reads 



u sp _ 

~V~ ~ 

where 



d 3 k 



In [2 cosh(/3£ fe /2)] 



hp(e Q - 2n)4 , 
(9) 

(10) 



are the Bogoliubov-likc elementary excitations with = 
h 2 k 2 /(2m) — /i. The saddle-point thermodynamic grand 
potential Q sp is related to the saddle-point partition func- 
tion Z sp and to the saddle-point effective action S sp by 
the formula 

Z sp = exp (-S ap /H) = exp (-/3 Q sp ) . (11) 

The value of </>o still needs to be determined. It is found 
from extremizing the action function S sp , or equivalcntly 
£l sp i with respect to <j) to determine the saddle-point 
value. Doing so one finds 



e - 2(i f d 3 k tanh (j3E k /2) 



(2tt) 3 2E k 



(12) 



The integral on the right side of Eq. (|12j) is formally 
divergent. However, expressing the bare detuning pa- 
rameter e in terms of the effective scattering length a s 
with the formula 1221 



e - 2(i 



d 3 k 



g 2 A<Kh 2 a s J (2Ttfh 2 k 2 

one obtains the regularized equation 

d 3 k /tanh(/3£ fc /2) m 



47rft 2 a, 



(2tt) 3 V 2£ fc 



(13) 



(14) 



The total number density n of atoms in the system is 
instead obtained from the thermodynamic formula n = 
—dfl sp /d(i, which gives 



d 3 k 

(2tt) 3 



[1 - (u 2 - vl) tanh ((3E k m] + 2$ . (15) 



with u 2 = (1 + £ fe /£fc) /2 and v 2 = (1 - $ k /E k ) /2. Eq. 



(|15p can be written as n = np + where ns = 2(f) 2 , is 
the density of atoms in the bosonic molecules (with two 
atoms per molecule). We stress that Eqs. (fT4|) and (fl"5|) 
generalize the gap and number equations one finds for a 
broad Feshbach resonance. In fact, setting 



A = gcj>o 



(16) 



the broad-resonance regime is easily obtained in the limit 
g — > +oo with the constraint of a finite energy gap A 
which implies 0o ~> 0. Notice that in the literature it 
is common to study models that in addition to the two- 
channel Feshbach resonant interaction considered above, 
a featureless nonresonant four-Fermi atomic interaction 
is also included [ll|. Gurarie and Radzihovsky [l|| have 
shown that in three dimensions doing so does not add any 
new physics to the pure two-channel model considered 
here. 



III. CONDENSATION 

The total condensate density no of the system is the 
sum of two contributions, namely no = riF,o+nB,Q, where 
Uf,o = 2|(V>|-(r,T)'i/;4_(r,T))| 2 is the condensate density 
of Cooper-paired atoms and Ub,o = 2\{<j>(r, r))| 2 is the 
condensate density of atoms in the preformed bosonic 
molecules. Within the saddle-point approximation we 
have n^o = 20q an d taking into account previous results 
on t ne condensate fraction of Cooper-paired atoms 
we finally obtain 



n 



d 3 k 

J2nf 



[u 2 k vl tanh 2 (pE k /2)] +2$. (17) 
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FIG. 1: (Color online). Upper panel: Scaled chemical poten- 
tial (J./tF as a function of the scaled inverse scattering length 
l/^Ffls). Lower panel: Scaled energy gap (iasa function of 
the scaled inverse scattering length l/{kFa s ). The curves cor- 
respond to different values of the scaled atom-molecule cou- 
pling 7 = ,gn 1//2 /eF- 



As discussed in [19j, at zero temperature the saddle- 
point approximation gives reliable results while at finite 
temperature it is necessary to include fluctuations about 
the saddle point cf>o: fluctuations are particularly impor- 
tant close to the critical temperature T c of the super- 
to-normal phase transition. For this reason we analyze 
Cooper pairing and Bose-Einstcin condensation in nar- 
row Feshbach resonances by using Eqs. (fll)) . (fT5j) and 
(fT7|) only at zero temperature, where they reduce to the 
following ones: 



V 
1 

no 
n 



3 r 



(18) 
(19) 



2V2 V A 



a A 2 

+ (20) 



where jl = [ijtp is the adimensional chemical poten- 
tial with cf = h 2 kp/(2m) the Fermi energy and kp = 
(37r 2 n) 1//3 is the Fermi wavenumber, A is the adimen- 
sional energy gap, y = \/{kpa s ) is the adimensional in- 
verse scatering length, 7 = gn 1 / 2 /ep is the adimensional 
resonant atom-molecule coupling, I\(x) and l2(x) are the 
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FIG. 2: (Color online). Upper panel: Total condensate 
fraction no/n as a function of the scaled inverse scattering 
length l/(fci?as). Middle panel: Condensate fraction n^o/n 
of Cooper paired atoms as a function of the scaled inverse 
scattering length l/(fc_Fa s ). Lower panel: Condensate fraction 
riBfi/n of atoms in preformed molecules as a function of the 
scaled inverse scattering length l/(fc_Fa s ). The curves corre- 
spond to different values of the scaled atom-molecule coupling 
7 = gn 1 



,V 2 



two monotonic functions 



h{x) 
h{x) 



+00 



+00 



y/(z 2 -x) 2 + l Z< 



^(z 2 - x) 2 + 1 



dz , (21) 



dz , (22) 



which can be expressed in terms of elliptic integrals. 

We have solved Eqs. (|T8|) and (jT9]) numerically obtain- 
ing the scaled (adimensional) chemical potential fi/eF 
and the scaled energy gap A/cf as a function of y = 
l/(kpa s ) for different values of 7 = gn 1 / 2 jtp. The re- 
sults are shown in Fig. [1] Remarkably the chemical 
potential [i does not depend very much on the resonant 
atom-molecule coupling 7 apart close to the unitarity re- 
gion (— 1 < y < 1). Instead the behavior the energy gap 
A depends strongly on 7 in the BEC region (y > 1). 
In particular, in the limit y — > +00 the energy gap A 
goes to a value which is obtained as the solution of Eq. 
(fTO]) with 7 2 (-oo) = 0, namely A/e F = j/V2- Taking 
into account Eq. (|16[) , this result means that in the limit 
y — > +00 one always has ng ~ n, for any finite value of 



7- 



In the upper panel of Fig. [5] we plot the total conden- 
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sate fraction no / n, obtained from Eq. ((20)) , as a function 
of the scaled inverse scattering length l/(kpa s ) for four 
values of 7. For a very small atom-molecule coupling, i.e. 
7 = 0.1 (dot-dot-dashed curves), the total condensate 
fraction becomes equal to 1 at unitarity y = 0. Indeed in 
this regime of weak-coupling the fraction no jn is mainly 
due to preformed molecules, as clearly shown in the mid- 
dle panel, where we plot the condensate fraction riF,o/n 
of Cooper paired atoms, and in the lower panel, where 
we plot the condensate fraction ns,o/n of atoms in pre- 
formed molecules. By increasing the resonant coupling 7, 
see for instance dot-dashed (7 = 1) and dashed (7 = 10) 
curves, the value of the total condensate fraction at uni- 
tarity (y = 0) decreases. Moreover, by increasing 7 in 
the BEC region (y > 1) it appears a finite condensate 
fraction of Cooper-paired (and bound) atoms while the 
fraction of condensed preformed molecules is strongly re- 
duced. Finally, for very large values of 7, e.g. 7 = 50 
(solid curves) of Fig. [3J in practice (i.e. in a large range 
of y) one recovers the broad-resonance limit where there 
are not preformed molecules and the BCS-BEC crossover 
is enterely due to Cooper-paired atoms. 

Within the zero-temperature saddle-point approxi- 
mation we have used, the condensate density riFfi 
of fermions is smaller than the total density uf of 
fermions, while the condensate density n^o of preformed 
molecules coincides with the total density ng of pre- 
formed molecules. In general, as shown by Werner, Tar- 
ruell, and Castin 24J , for the two-channel model the den- 



sity riB (but not ub,o nor uf,o) can be directly related to 
the derivative of the energy of the gas with respect to the 
inverse scattering length. The zero-temperature conden- 
sate depletion (jib — riB,o)/ n B of bosonic molecules can 
be instead obtained going beyond the saddle-point ap- 



proximation by including quadratic fluctuations of <j)(r, t) 
around O in the action functional. 



IV. CONCLUSIONS 

In this paper we have used the two-channel model of 
atomic fermions coupled to preformed bosonic molecules 
to study two crossovers: the BCS-BEC crossover driven 
by the s-wave scattering length and the crossover from 
a narrow to a broad resonance controlled by the atom- 
molecule coupling. Within the saddle-point approxima- 
tion of the path integral formalism we have derived an 
extended gap equation which depends on both the s-wave 
scattering length and the atom-molecule coupling. Only 
for a very large atom-molecule coupling the extended 
gap equation reduces to the one of broad Feshbach res- 
onances. At zero temperature, we have calculated the 
chemical potential and the energy gap of the system. 
Remarkably, in the deep BEC region of the BCS-BEC 
crossover the energy gap A goes to the asymptotic value 
jn 1 / 2 /^, where g the atom-molecule coupling and n 
is the total atomic density. We have then analyzed the 
condensate density and the condensate fraction of the 
system. We have found that by increasing the resonant 
atom-molecule coupling the condensate fraction of pre- 
formed molecules is strongly reduced in the BEC regime 
with a corresponding growth in the condensate fraction 
of Cooper-paired atoms. As expected, for very large 
values of the resonant coupling one recovers the broad- 
resonance regime characterized by the absence of pre- 
formed molecules. 
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